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Abstract. Recent studies have claimed the existence of very massive stars (VMS) up to 300 A/0 
in the local Universe. As this finding may represent a paradigm shift for the canonical stellar 
upper-mass limit of 150 Mq, it is timely to discuss the status of the data, as well as the far- 
reaching implications of such objects. We held a Joint Discussion at the General Assembly in 
Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the 
formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim 
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was to reach broad consensus between observers and theorists on how to identify and quantify 
the dominant physical processes. 

Keywords. Keywordl, keyword2, keyword3, etc. 



1. Introduction 

The last decade has seen a growing interest in the study of the most massive stars, as 
their formation seems to be favourable in the early Universe at low metalliciy (Z), and is 
thought to involve a population of objects in the range 100-300 M© (Bromm et al. 1999; 
Abel et al. 2002). The first couple of stellar generations may be good candidates for the 
reionization of the Universe. Notwithstanding the role of the first stars, the interest in 
the current generation of massive stars has grown as well. Massive stars are important 
drivers for the evolution of galaxies, as the prime contributors to the chemical and energy 
input into the interstellar medium (ISM) through stellar winds and supernovae (SNe). 
A number of exciting developments have recently taken place, such as the detection of a 
long-duration gamma-ray burst (GRB) at a redshift of 9.4, just a few hundred millions 
years after the Big Bang (Cucchiara et al. 2011). This provides convincing evidence that 
massive stars are able to form and die massive when the Universe was not yet enriched. 

The specific reason for holding this JD was the recent evidence for the existence and 
subsequent deaths of very massive stars (VMS) up to 300 M©. Gal- Yam et al. (2009) 
claimed the detection of a pair-instability SN (PSN) from a VMS. These explosions 
are thought to disrupt stars without leaving any remnants. Crowther et al. (2010) re- 
analyzed the most massive hydrogen-and nitrogen-rich Wolf-Rayet (WNh) stars in the 
center of R136, the ionizing cluster of the Tarantula nebula in the Large Magellanic Cloud 
(LMC). The conclusion from their analysis was that stars usually assumed to be below the 
canonical stellar upper-mass limit of 150 M Q (of e.g. Figcr 2005), were actually found to 
be much more luminous, and with initial masses up to ^320 M©. Prior to discussing the 
formation, evolution, and fate of VMS, and before we should explore the full implications 
of these findings, it is imperative to discuss the various lines of evidence for and against 
VMS. 

VMS are usually found in and around young massive clusters, such as the Arches cluster 
in the Galactic centre and the local starburst region R136 in the LMC. Such clusters may 
harbor intermediate-mass black holes (IMBHs) with masses in the range of several 100 to 
several 1000 M© and may provide insight into the formation of supermassive black holes 
of order 10 5 M©. Young clusters are also relevant for the unsolved problem of massive 
star formation. 

For decades it was a struggle to form stars over 10 M©, as radiation pressure on dust 
grains might halt and reverse the accretion flow onto the central object (e.g., Yorke & 
Krucgcl 1977). Because of this problem, astrophysicists have been creative in forming 
massive stars via competitive accretion and merging in dense cluster environments (e.g., 
Bonncll et al. 1998). In more recent times several multi-D simulations have shown that 
massive stars might form via disk accretion after all (e.g., Krumholz et al. 2009; Kuiper 
et al. 2010). In the light of recent claims for the existence of VMS in dense clusters, 
however, we should redress the issue of forming VMS in extreme environments. 

Massive clusters may be so dense that their early evolution is largely affected by stellar 
dynamics, and possibly form very massive objects via runaway collisions (e.g., Portegies 
Zwart et al. 1999; Giirkan et al. 2004), leading to the formation of VMS up to 1000 M© 
at the cluster center, which may produce IMBHs at the end of their lives, but only if 
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VMS mass loss is not too severe (see Belkus et al. 2007, Yungelson et al. 2008, Glebbeek 
et al. 2009, Pauldrach et al. 2012). 

VMS are thought to evolve almost chemically homogeneously (Grafener et al. 2011), 
implying that knowing the exact details of the mixing processes (e.g., rotation, magnetic 
fields) could be less relevant in comparison to their canonical ^10-60 Mq counterparts. 
Instead, the evolution and death of VMS is likely dominated by mass loss. A crucial 
issue regards the relevance of episodes of super-Eddington, continuum-driven mass loss 
(such as occurs in Eta Carinae and other Luminous Blue Variable (LBV) star eruptions), 
which might be able to remove large amounts of mass - even in the absence of substantial 
line-driven winds (see Sect. [6]). 

A final issue concerns the fate of VMS, and more specifically whether VMS end their 
lives as canonical Wolf-Rayet (WR) stars giving rise to Type Ibc SNe, or do they explode 
prematurely during the LBV phase? Might some of the most massive stars even produce 
PSNe? And how do PSNe compare to the general population of super-luminous SNe 
(SLSNe) that have recently been unveiled by Quimby et al. (2011), and are now seen out 
to high rcdshifts (Cooke et al. 2012). Such spectacular events can only be understood 
once we have obtained a basic knowledge of the physics of VMS. 

In Section [2] a definition for VMS is adopted, before we discuss the evidence for and 
against super-canonical stars, i.e. objects above the traditional 150 Mq stellar upper 
mass limit (Sect. [3]). After it is concluded that VMS probably exist, the next question to 
address is how to form such objects (Sect. 2]). We then discuss the properties of VMS in 
Sect. 03 before discussing their mass loss (Sect. [HI, evolution and fate (Sect. [7]). We end 
with the implications (Sect. [8]) and final words. 



2. Definition of VMS 

Before we can discuss the evidences for and against VMS, one of the very first issues 
we discussed during the JD was what actually constitutes a "very" massive star. One 
may approach this in several different ways. Theoretically, "normal" massive stars with 
masses above ~8M Q are those that produce core-collapse SNe (Smartt et al. 2009), but 
what happens at the upper-mass end (Nomoto 2012)? 

Above a certain critical mass, one would expect the occurrence of PSNe, and ideally 
this could be the lower-mass limit for the definition of our VMS. However, in practice this 
number is not known a priori due to mass loss, and as a result the initial and final masses 
are not the same. In other words, the initial main-sequence mass for PSN formation is 
model-dependent, and thus somewhat arbitrary. Moreover, there is the complicating issue 
of pulsational pair-instability supernova (Puls-PSN) at masses below those of full-fcdgcd 
PSNe (e.g. Woosley et al. 2007). One could of course resort to the mass of the helium (He) 
core for which objects reach the conditions of electron/positron pair-formation instability. 
Heger showed this minimum mass to be ^40 Mq to encounter Puls-PSN and ~65 Mq to 
encounter PSNe (see also Chatzopoulus & Wheeler 2012). 

An alternative definition could involve the spectroscopic transition between normal 
main-sequence O-type stars and hydrogen- rich Wolf-Rayet stars (of WNh type), which 
have also been shown to be core H burning main-sequence objects. However, such a 
definition would also be model dependent. 

For these reasons, we took the decision to follow a more pragmatic approach: we 
consider stars to be very massive when the initial masses are ^100 Mq, or higher. 
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Figure 1. VLT MAD K s -band 12 x 12 arcsec (3x3 parsec for the LMC distance of 49 kpc) 
image of R136 (Campbell et al. (2010) in conjunction with a view of the central 4x4 arcsec (1 
x 1 parsec) in which the very massive WN5h stars discussed are labelled (component b is a lower 
mass WN9h star). Relative photometry agrees closely with integral field S1NFON1 observations 
of Schnurr et al. 2009). See Crowther et al. (2010) for details. 



3. On the existence of VMS 

With this definition, the question of whether very massive stars exist can convincingly 
be addressed with the one answer: yes they do! However, the pertinent issue for our JD 
was whether the widely held "canonical" upper-mass limit of 150 Mq has recently been 
superseded. 

Paul Crowther gave the first invited review presenting evidence for initial masses as 
high as 320 Mq. Crowther provided historical context to some of the astronomical com- 
munity's sceptism regarding such high masses in R136, in particular the spectacular claim 
for the existence of a 2500 Mq star R136 in the 30 Doradus region of the LMC some 3 
decades ago (e.g. Cassinelli et al. 1981). Higher spatial resolution showed that R136 was 
actually not a single supermassive star, but the object broke up into a young cluster 
containing significantly lower mass objects, including the current record holder R136al 
(as well as R136a2, R136a3, and R136c). Therefore, over the last few decades, there was 
a consensus of a 150 Mq stellar upper mass limit (Weidner & Kroupa 2004; Figer 2005; 
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Oey & Clarke 2005, Koen 2006), although the accuracy of this magic number of 150 was 
low (Massey 2011). 

Crowther et al. (2010) re-analyzed the photometric and spectroscopic data of the VMS 
in R136. In comparison to the older WFPC2 data, they used ground-based adaptive optics 
photometry (see Fig.[l|. In combination with their spectral analysis using the CMFGEN 
non-LTE atmosphere code of Hillier & Miller (1998), this lead to higher estimates for 
effective temperatures and bolometric corrections. In conclusion they claimed that the 
R136 cluster hosts several stars with masses as high as 200-300 Mq. 

Crowther et al. also performed a "sanity check" on similar WNh objects in the Galactic 
starburst cluster NGC3601. Although these objects were fainter, and less massive, than 
those in R136, the advantage was the available dynamical mass estimate by Schnurr et al. 
(2008) of the binary object NGC3601-A1 of 116 ± 31 + 89 ± 16 Mq. This is important 
as the least model-dependent way to obtain stellar masses is through the analysis of the 
light-curves and radial velocities induced by binary motions. We also note that Rauw et 
al. (2004) and Bonanos et al. (2004) found both components of the eclipsing Wolf-Rayet 
binary WR20a to be particularly massive, with 83 ± 5 and 82 ± 5 Mq, with small error 
bars. 

During the lively discussion that followed Paul's review, some attendants argued that 
the luminosities derived by Crowther et al. are uncertain and these "single" objects might 
actually contain multiple sources. Whilst short-period binaries were not detected by 
Schnurr et al. (2009), longer period binaries are harder to exclude. One of the additional 
arguments by Crowther et al. was that X-rays have not been detected, whilst they may 
have been expected on the basis of colliding wind binary (CWB) simulations by Pittard & 
Stevens (2002). Oskinova however countered this argument on the basis that empirically 
a low X-ray luminosity cannot serve as a robust argument against a binary nature (see 
Sect.!]). 

Najarro noted that even the best image so far of the Arches cluster with Keck (Fig. [5]). 
has a limited spatial resolution of 50 milli-arcsec (mas), which corresponds to roughly 
1/10 of the diameter of the circles marking the PSF reference stars in the figure. Given 
that the LMC is almost 7 times further away, and that the VLT is smaller than Keck, the 
circles from Fig. [2] would roughly correspond to the spatial resolution achieved with the 
VLT if the Arches clusters was in the LMC. In other words, the Arches stars would effec- 
tively "merge" with surrounding objects. This analogy suggests that we cannot exclude 
the possibility that the bright WNh stars in the R136 cluster core could still break up 
into several lower-mass WNh stars. However, for example a ^300M Q star could at best 
break up into a pair of ~150M Q stars (given the shallow slope of the stellar luminosity 
to mass ratio at the high-mass end). 

Moreover, Bestenlehner noted that there is a near-identical twin of R136a3 WNh star 
in 30 Doradus: VFTS 682. Its key relevance is that it is found in isolation from the R136 
cluster, some 30 pc away. For this object line-of-sight contamination is far less likely 
than for the R136 core stars. VFTS 682 thus offers another sanity check on the reliability 
of the luminosities for the R136 core stars. Bestenlehner et al. (2011) derived a high 
luminosity of log(L/L Q ) = 6.5 and a present-day mass of 150 Mq for VFTS 682. This 
implied an initial mass on the zero-age main sequence (ZAMS) higher than the canonical 
upper-mass limit - likely ~ 200 Mq . 

In his talk Hamann showed results from another extremely luminous WN star in the 
Galactic center region: the Peony star (Barniske et al. 2008). The luminostiy of this star 
is determined from spectral analysis as log(L/LQ) = 6.5 ± 0.2 and initial mass between 
150 ^ M u /Mq ^ 200. The star is located above the Humphreys-Davidson limit, in the 
region populated by the LBV stars. However, the hydrogen content is lower in WR 102ka 
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Figure 2. K-band mosaic of the core field of the Galactic Arches cluster at a resolution of 50 
mas and with a sensitivity of Kn m = 20.6mag. The figure is taken from Clarkson et al. (2012). 

compared to the Pistol star, while helium is higher. This indicates a more advanced 
evolutionary stage of the former, compared to "normal" LBV stars. 

In summary, whilst one cannot exclude that the object R136al claimed to be ^300 M Q 
in the R136 cluster might still turn out to "dissolve" when higher spatial resolution ob- 
servations were to become available, a number of sanity checks involving binary dynamics 
and "isolated" objects make it rather convincing that stars with ZAMS masses of 200 M Q 
exist. In any case, no fundamental reason was identified why "150" would be a magic 
number. 



4. Formation of VMS 

The key question to address was what is so special about the formation of 250 M stars 
in comparison to 'normal' 50 M Q O-type stars? Can these VMS only form inside dense 
cluster environments? Other relevant questions involve the time-scale and evolutionary 
stages when they finish formation/accretion. How fast do they rotate? What VMS frac- 
tion is in binaries, and what are the binary properties? Are the most massive stars all 
the result of stellar mergers? 
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4.1. Theoretical considerations 



According to our second invited review speaker, Mark Krumholz, VMS formation is 
not a fundamentally different problem from the formation of massive stars in general. 
From an observational standpoint VMS in clusters appear as part of a continuous initial 
mass function (IMF), with no special features that mark them off as different from the 
remainder of the stellar population. From a theoretical standpoint, both the radiative 
and wind luminosities of stars are increasing functions of mass, but from the standpoint 
of star formation there is no natural dividing line that would put a star of, for example, 
50 Mq into one category and a star of 250 Mq into another. For this reason, it makes 
sense not to treat the formation of very massive stars as a separate problem, and instead 
to embed it in the broader context of forming the upper end of the IMF. 

There are three main challenges to forming massive stars: fragmentation, binarity, 
and radiation pressure. The fragmentation problem is simple to state: we certainly see 
gas over-densities ("cores") with masses of ~ 100 Mq or larger a nd radii of ~ 0.1 pc 
or less that seem pro mising sites for massive star formation (e.g. iBeuther et al. 2005 : 
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2010). However, given these masses and radii, and typical molecular 
cloud temperatures of ~ 10 K, the Jeans mass is below 1 Mq, so how can the objects 
avoid fragmenting and collapse to form an object containing 100 Jeans masses or more? 
The answer seems to be that the classical Jeans analysis of an isothermal gas does a rather 
poor job at predicting the behavior of a radiatively heated, magnetized fluid. Once small 
stars form and begin to accrete, their radiation heats the gas around them and suppresse s 
fragmentation in it (jKrumholzj 120061: iKrumholz et all 120071: IKrumholz fc McKeel l2008h 
Magn etic fields also make it much more difficult for the gas to fragment ( Hennebelle et al.l 
20111 ). and the combined e ffects of radiation and magnetic fields seems to be particularly 
effective, as seen in Fig.l3l ( Commercon et al.l 1201 It iMvers et al.ll2012l ). These effects to- 
gether seem to resolve the fragmentation problem, indicating that massive stars can form 
from the direct collapse of massive protostellar cores with properties similar to those ob- 
served. 

The second challenge is explaining why so many massive stars appear to be binaries 
(e.g. Sana et al. 2012). The observed multiplicity fraction rises very sharp ly with stella r 



mass, reaching near unity for O and B stars (excluding runaways - e.g. lBrownll200lf ) 



Radiation-hydrodynamic simulations of star cluster format ion appear able to replicate 



the o bserved dependence of multiplicity on primary mass (jBatd 12012c IKrumholz et al 



2012h . and digging into the physical origin for this result indicates that it combines 
two effects. The first is simple N-body processing: close encounters between stars in 
young clusters tend to put the most massive members into binaries even if they are not 
born that way, while strippin g companions from less m assi ve stars. The second is dis k 
fragmentation, as explored by Kratter fc Matzner ( 200(jh and Kratter et al ] (|2008ll2010h . 
Massive stars form with high accretion rates, and these accretion rates tend to produce 
disks with masses that can approach that of the primary. When this happens, disks are 
likely to fragment, and a common outcome of this process is that the disk produces a 
massive companion to the primary. 

The third challenge is radiation pressure. Dusty interstellar gas has a high opac- 
ity, and the radiation force is proportional to this opacity. As a result, the radiative 
force exert ed by a star's light can exceed its gravitational fo rce for all stars larger than 
- 20 M Q (|Wolfire fc Cassinellil Il987t IKrumholz etaH l2009h . However, recent numeri- 
cal and analytic work has shown that this problem is mostly an illusion; the radia- 
tion pressure barrier can be circumvented in numerous ways. First, in the presence of 
an accretion disk, radiation can be beamed away from the bulk of the accreting mat- 
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Figur e 3. Column den sity projections from three simulations of a collapsing high mass core, 
from lMvers et all (|2012T ). All three simulations use identical initial conditions, but the leftmost 
one includes radiation and magnetic fields, the center one includes magnetic fields but not 
radiative transfer, and the rightmost one includes radiation but not magnetic fields. Note the 
dramatic reduction in number of stars (black circles) in the radiation plus magnetic fields case. 
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ter (iNakano et alJ ll995l IJiiina fc Adamd 119961 lYorke fc Sonnhalterl [20021 iKuiper et al 



20101 120111 ). Second, radiation-driven Ray leigh- Taylor instabilities can break up bubbles 



of radiation and allow matter to a ccrete through optically-thick fingers (jKrumholz et al 



2009t IJacquet fc KrrjmholjlioTlh . Figure H shows an example. Third, protostellar jets 
can punch holes in accreting cores that allow radiation to leak out, reducing the net 
radiation force in most directions and allowing accretion to con tinue even though, av- 
eraged over 47T sr, the r adiation force is stronger than gravity ( Krumholz et al.l 12005 ; 



Cunningham et al.l 1201 If) . The takeaway message from all of these simulations is that 
radiation pressure poses no barrier to the formation of stars up to arbitrary masses. 

Given that modern theoretical models have removed all of the serious objections to 
forming massive stars in the same way that low mass stars form, i.e. by accretion from a 
collapsing gas core, there is no need to resort to exotic processes like stellar collisions (e.g. 
Bonnell et al.lll998f ). This does not mean, however, that stellar collisions cannot happen 
and also contribute to the massive star population. The question is: under what circum- 
stances do we expect collisions to be important? The mos t rec ent and comprehensive 
paper s to address this question are iMoeckel fc Clarke (|201ll ) and iBaumgardt fc Klessen 
(|2011l ) who both conducted N-body simulations of stars confined by a gaseous potential, 
and who came to similar conclusions. They find that the collisional formation of very 
massive stars is only significant if the stellar density is extremely high, in excess of 10 7 
pc -3 , with surface densities reaching 10 5 pc -2 . These numbers are so high that even the 
Arches cluster would not be expected to have significant contributions to its massive star 
population by collisions. Moreover, if collisions are significant, they do not produce an 
IMF that looks like what we observe. Instead, because collisions tend to occur among the 
most massive stars, they produce an IMF that has a Salpeter-like slope at intermediate 
masses, then a deficit of stars at higher masses, and finally one or a few very massive 
stars. IMFs with dips of this sort have not been observed, again suggesting that collisions 
are not likely to be significant contributors to the massive star populations. 

4.2. The potential role of dynamically induced mergers 

An alternative VMS formation scenario was presented by Sambaran Banerjee (with Pavel 
Kroupa and Seungkyung Oh; Banerjee et al. 2012a) who argued that super-canonical 
stars can be formed out of a dense stellar population - with a canonical IMF and with a 
150 Mq upper limit - through dynamically induced mergers of the most massive binaries. 
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Figure 4. Volume renderings of a simulation from lKrumholz et al.l (|2009T ) involving the density 

field in a (4000AU) 3 region at 55.0 kyr of evolution. The color scale is logarithmic and runs from 
10 -i6.5 _ 10 -i4 

g cm . The left panel shows the polar view, and the right one denotes the 
edge-on view. The figure highlights how Rayleigh- Taylor instability fingers channel matter onto 
a massive binary star system. 

Banerjee et al. performed direct N-body computations (NBODY6; Aarseth 2003) of a 
fully mass-segregated star cluster mimicking R136 in which all the massive stars are in 
primordial binaries. Banerjee et al. account for the mass evolution of the super-canonical 
stars and the resulting shortened (« 1.5 Myr) lifetimes in their super-canonical phases 
using stellar evolutionary models by Kohler & Langer (2012) that incorporate Vink et 
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al. (2000) mass-loss rates for the main sequence and Hamann et al. (1995) for the He 
burning WR phase. 

Banerjee et al. find that super-canonical stars begin to form via dynamical mergers 
of massive binaries from « 1 Myr cluster age, obtaining stars with initial masses up 
to w 250 Mq. Multiple super-canonical stars are found to remain bound to the cluster 
simultaneously within a super-canonical lifetime. Banerjee also noted that some of these 
objects can be formed at runaway velocities which escape the cluster at birth. For in- 
stance, the most massive apparently isolated WNh star VFTS 682 might be an expected 
slow runaway (Banerjee et al. 2012b; see also Fujii & Portegies-Zwart 2012). 

The Banarjee et al. models indicate that had super-canonical stars formed primordially 
alongside the rest of the R136 cluster, i.e. violating the canonical upper limit, they would 
have evolved below the canonical 150 Mq limit by « 3 Myr, the likely age of the bulk 
of R136 according to Andersen et al. (2009). In other words, Banerjee et al. argue that 
primordially-formed super-canonical stars should not be observed at the present time in 
R136, whilst it is quite plausible that a collection of dynamically formed super-canonical 
VMS would be observed in the centres of young massive starburst clusters. 

A fully self-consistent N-body computation incorporating detailed accurate evolution- 
ary and mass-loss recipes would be needed to confirm these scenarios. 

4.3. Rotation rates as a constraint on massive star formation 

Returning to the more conventional ways of forming massive stars, Anna Rosen addressed 
the question of what sets the initial rotation rates of massive stars. The physical mech- 
anisms that set the initial massive star rotation rates are a crucial unknown in current 
star-formation theory. Observations of young, massive stars provide evidence that they 
form in a similar fashion to their lower mass counterparts. The magnetic coupling between 
a star and its accretion disk may be sufficient to spin down low-mass pre-main-sequence 
(PMS) stars to well below breakup at the end stage of their formation when the accretion 
rate is low. However, Anna showed that these magnetic torques are insufficient to spin 
down massive PMS stars due to their short formation times and high accretion rates. 
Anna developed a model for the angular momentum evolution of stars over a wide range 
in mass, considering both magnetic and gravitational torques. She finds that magnetic 
torques are unable to spin down either low-mass or high-mass stars during the main 
accretion phase, and that massive stars cannot be spun down significantly by magnetic 
torques during the end stage of their formation either. Spin-down occurs only if massive 
stars' disk lifetimes are substantially longer or their magnetic fields are much stronger 
than current observations suggest (Rosen et al. 2012). 

4.4. Observations of massive star formation 

Heavy extinction hides the birthplaces of massive stars from view, and the short formation 
timescales set a strong limitation to the sample of objects that can be studied. So far, 
our physical knowledge of massive young stellar objects (YSOs) has been derived from 
near-IR imaging and spectroscopy, revealing populations of young OB-type stars, some 
still surrounded by a disk, others apparently 'normal' main sequence stars powering H II 
regions. The most important spectral features of OB-type stars are however located in 
the ultraviolet (UV) and optical range. 

Kaper (with Ellerbroek, Ochsendorf and Bik) showed that with the new optical/near- 
infrared spectrograph X-shooter on ESO's Very Large Telescope (VLT), it is possible to 
extend the spectral coverage of these massive YSOs into the optical range. First results 
are very promising, although they seem to probe the intermediate- mass (^2 — 8 M ) 
range rather than the massive star range. Ellerbroek et al. (2011) discovered a jet (HH 
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Figure 5. The location of B275 (red parallelogram) in the HRD next to PMS tracks from 
Hosokawa et al. (2010) with the ZAMS mass labeled and open symbols indicating lifetimes. The 
thin dashed and thin dot-dashed lines are the birth lines for accretion rates of 10 -4 Mq yr _1 
and 10" 5 M yr" 1 , respectively; the thick solid line is the ZAMS (Schaller et al. 1992). The 
filled and open circles represent stars in M 17 for which a spectral type has been determined 
(Hoffmeister et al. 2008); dots are other stars in M 17. B275 is on its way to becoming a 6 8 Mq 
ZAMS star, so far one of the most massive pre-main-sequence star known. The figure has been 
adapted from Ochsendorf et al. (2011). 
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1042) produced by the massive YSO nr292 in the massive star forming region RCW36, 
demonstrating that the object is still actively accreting. The mass of this star is likely less 
than 6 M Q and remains uncertain as photospheric features are not detected. The first 
firm spectral classification of B275, a massive YSO in M17, results in its precise location 
on a PMS track for a ~ 7M Q star (see Fig.[5J Ochsendorf et al. 2011), has the size of a 
bloated giant, as predicted by models of Hosokawa et al. (2010), and is still surrounded 
by a disk. 

It remains unclear whether the progenitors of the most massive VMS will be detectable 
this way. 

4.4.1. Radiation Transfer Modeling: From massive to very massive star formation 

In order to properly interpret the observations of massive protostars, Zhang et al. (with 
Tan, McKee, and de Buizcr) presented a radiation transfer (RT) model for a massive core 
in a high pressure environment that forms a massive star through core accretion (Zhang 
& Tan 2011; Zhang, Tan & McKee, 2012, ApJ submitted). This RT model is based on 
the Turbulent Core model of McKee & Tan (2003): massive stars form in dense clumps 
with high surface density. Assuming the rotational-to-gravitational energy ratio is 2% in 
the core, whilst the disk has a diameter of ^1000 AU. The high accretion rate leads to 
a disk mass comparable to the stellar mass. In Zhang et al. (2012), the treatment of the 
disk was improved by allowing radially varying accretion rates due to a supply of mass 
and angular momentum from the infall envelope and their loss to the disk wind. The 
transfer of accretion power to mechanical power of the wind was also accounted for. An 
approximate disk wind solution was developed partly based on the Blandford & Payne 
(1982) model. The simulation was performed with the latest version of the Monte Carlo 
RT code by Whitney et al. (2003). Corrections made by adiabatic cooling/heating and 
advection were included, and so were the gas opacities. 

The model was compared to the massive protostar G35.2-0.74N. At a distance of 2.2 
kpc, radio continuum emission indicates there is a bipolar outflow from this source. In 
recent SOFIA-FORCAST observations of the massive protostar G35.2-0.74N at 31 and 
37m (Zhang et al. 2012, in prep.) both the near- and far-facing sides of the outflow can 
be seen. The latter is missing in mid-IR continuum. By fitting both the observed SED 
and the outflow-axis intensity profiles with the RT model, the bolometric luminosity was 
inferred to be ~ 10 Lq after correcting for foreground extinction and the dependence of 
luminosity on inclination (flashlight effect). 

The fitting model suggests a massive ~ 30 M Q protostar is forming from a ~240 M Q 
core with a high surface density of ~lg/cm 2 , via relatively ordered collapse and accretion, 
and driving powerful bipolar outflows. These results seem to support the core accretion 
theory which predicts massive stars may form in similar ways to their low-mass coun- 
terparts. Simply extending the model to protostars with higher masses up to 100 Mq, 
Zhang et al. noticed a shift of the SED peaks to shorter wavelengths and a change of 
mid-IR slope. The flashlight effect turns out to be huge for such sources. In the mid-IR, 
the luminosity from face-on view can be higher than that from edge-on view by as much 
as 3 dex - suggesting that a large correction factor needs to be applied when using IR 
observations for deducing source luminosities. 

5. Properties of VMS 

Once a star has formed it starts to burn hydrogen on the ZAMS. For massive stars, the 
main sequence band probably consists of two spectroscopically distinct groups of objects: 
the O stars, and the WNh stars. The first group are thought to form the normal massive 
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star sequence with masses of up to at least ^6OAf . For the most massive "O stars on 
steroids" WNh (H-rich WR stars), their masses may be up to 300 M . or higher. 

The transition mass between O and WNh is not exactly known but is probably lies 
somewhere in the range 60-120AfQ. Of course these numbers are model (mass loss and 
Z) dependent. The WNh stars are presumably in close proximity to the Eddington limit, 
and so are their descendant LBV and classical WN, WC, WO Wolf-Rayet stars. Some of 
these objects are thought to be in an evolutionary phase just prior to their final demise. 

At the end of this section, we also consider extra-galactic VMS properties, both at 
very low (e.g. IC 1613) and high metal content (e.g. M33). 

5.1. WNh H-rich Wolf-Rayet stars - "0 stars on steroids" 

WNh stars were discussed by Crowther and Hamann. In their analysis they make use 
of non-LTE model atmospheres, such as CMFGEN (Hillier & Miller 1998) and PoWR 
(Hamann et al. 2006), whilst FASTWIND (Puis et al. 2005) is often used for normal 
O stars with weaker winds. These are ID spherical symmetric codes that include line- 
blanketing and stellar winds allowing for micro-clumping (optically thick macro-clumping 
allowing for porosity is non-standard). Another important aspect of the analysis is that 
of accurate infrared photometry (see the infrared SpS on atmosphere modelling). This is 
especially true for those objects in the cores of dense clusters such as R136. The upshot of 
recent analyses is that WNh stars are more luminous, and more massive, than previously 
thought (see Sect. [3]). Hamann presented some early results from their spectral analysis 
of WN stars in the Magellanic Clouds (Hainich et al. in prep.) - complementing the 
published Galactic study of Hamann et al. (2006). Whilst the LMC objects are found 
to be bright with logL/ L@ values in the range 5.3-5.9, Hamann warns that many of the 
brighter allegedly single WNh stars may actually be binaries, and not accounting for this 
fact may overestimate the luminosities of WNh stars. 

5.2. LBVs: unstable massive stars close to the Eddington limit 

Very massive stars are thought to evolve through the unstable Luminous Blue Vari- 
able phase, when enormous amounts of mass are lost. While LBVs have been classically 
thought to be rapidly evolving massive stars in the transitory phase from O-type to Wolf- 
Rayet stars, recent studies have suggested that LBVs might surprisingly explode prema- 
turely as a core-collapse supernova (Kotak & Vink 2006, Gal- Yam et al. 2007, Mauerhan 
et al. 2012). Such a striking result highlights that the evolution of VMS through the LBV 
phase is far from being understood. 

Groh discussed the recent advances in understanding LBVs, in particular how to dis- 
tinguish them from the normal B supergiants and hypergiants (as discussed by Negueru- 
ela). LBVs can be recognized either from Giant Eruptions like Eta Car and P Cygni, or 
through their S-Dor variability, sometimes imprinted on peculiar looking double-peaked 
absorption profiles, in just a single epoch spectrum (Groh & Vink 2011). LBVs do not 
always present both S Dor and giant outburst phenomena, leaving room for quite a het- 
erogeneous class of objects. 

Groh emphasized that LBVs are characterized on a phenomenological basis and, there- 
fore, LBVs are neither a spectroscopic nor an evolutionary classification. Particular em- 
phasis was given to describe the main properties of the S-Dor type variability (also 
performed by Stringfellow). These are changes in the hydrostatic radius and bolometric 
luminosity. Finally, the role of rapid rotation on LBVs was discussed (Groh et al. 2009). 
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5.3. The Galactic WC stars 

WC star spectra are well known for their broad emission lines from helium, carbon and 
oxygen. Due to the absence of hydrogen these stars have to be core-helium burning and 
have as such been identified as the late evolutionary stage of massive and very massive 
stars. Sander et al. (2012) analyzed the optical and UV spectra of over 50 Galactic WC 
single stars and derived their parameters including their mass-loss rates. Sander et al. 
showed that the positions of the Galactic WC stars in the HR-diagram do not fit with 
the assumption that the most massive stars will pass the WC stage. Instead Sander et al. 
argue that their results indicate that WC stars come from an initial mass range between 
20 and somewhere around 50 Mq. The stellar evolution models of Vanbeveren et al. 
(1998) that include enhanced mass loss during the red sugergiant (RSG) phase appear 
to properly account for the location of the WC stars in the observed HR-diagram. It 
also seems that stars with higher initial masses do not reach the WC stage but instead 
explode after passing a WNL and probably an LBV stage. 

5.4. WO stars 

Sander et al. also showed that the Galactic W02 stars WR 142 and WR 102 have signif- 
icantly different parameters from the WC stars. The WO positions in the HR diagram 
suggest that these stars could be close to or already in the stage of carbon burning which 
makes them interesting SN Ic candidates (Georgy et al. 2012; Yoon et al. 2012). 

Tramper et al. presented VLT/X-Shooter spectroscopy of DR1, a W03 star in the low- 
metallicity galaxy IC 1613. A preliminary spectroscopic analyis using CMFGEN indicates 
a high temperature (of ~150 kK) which is also supported by the very strong nebular He 
II emission. The oxygen abundance does not seem to be enhanced compared to values 
found for early-type WC stars, suggesting that the strong oxygen emission is likely a 
temperature effect, rather than being caused by an increased oxygen abundance. 

5.5. Rapidly rotating WR stars 

The WR rotation issue is especially relevant in view of the suggested link between rotating 
WR stars and long-duration GRBs. WR 2 (WN2) is a well known oddball. It is the most 
compact and the hottest (~ 140 kK) pop I WN star known in the Galaxy. It is also 
one of the best candidate for strange- mode pulsations (Glatzel et al. 1999), but they 
are not observed in photometric observations. The spectrum displays bowler-hat shaped 
emission lines, in contrast with the more normal Gaussian or even flat-top and triangular 
profiles of most other WR winds. Remarkably, Hamann et al. (2006) have analyzed WR 

2 using their latest model-atmosphere code, and the model spectrum failed to reproduce 
its line-profiles, unless it was folded with a rotation curve near the break-up limit, i.e. 
1900 km/s. 

However, Chene presented the polarized spectrum of WR 2, which shows no sign of 
wind asymmetry expected for such rapid rotation. Interestingly, WR 2 appears to dis- 
play clumps that are moving in a similar fashion as in the optically thinner wind of WR 

3 (Chene et al. 2008). Hence, the shape of WR 2 spectra line cannot be the result of 
extreme opacity either. 



5.6. Massive stars in the low metallicity galaxy IC1613 

Because low metallicity environments may favour the formation of massive and very 
massive stars, we switch our attention to extra-galactic properties. Herrero discussed the 
selection of massive star candidate stars in the low-Z galaxy IC 1613 from the catalog 
of Garcia et al (2009). IC 1613 has a metallicity Z= 0.13 Z© according to the analysis 
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of B-supergiants by Bresolin et al. (2007), or Z= 0.08—0.15 Z Q according to several 
analyses of HII regions (see Herrero et al. 2012). Herrero et al. observed the selected 
stars with OSIRIS/GTC at a resolution R = 1000 to determine their spectral types, as 
a previous step for a more detailed analysis. They were able to classify 12 new OB stars 
and confirmed one more known O-type star. 

The spectra of the O stars were good enough for quantitative analysis (albeit with 
errors slightly larger than in typical analyses, see Repolust et al. 2004). Combining these 
results to those from the literature (Tramper et al. 2011; Herrero et al. 2012), Herrero 
presented the first effective temperature scale for sub-SMC metallicities. This tempera- 
ture scale is slightly hotter than that derived for SMC stars from the data of Massey et 
al. (2009), Mokiem et al. (2007), and Trundle et al. (2007). 

5.7. Stellar abundances from massive stars in M33 

With its distance less than a Mpc (Bonanos et al. 2006) and its favourable inclination 
angle, this makes M33 an ideal galaxy to study the chemical evolution in spiral galaxies. 
Not long ago, the only way to carry out detailed chemical analyses of nearby galaxies 
was through the quantitative studies of H II regions. However, the abundances derived 
from massive OB stars are the tracers of the present-day chemical composition, providing 
information that cannot be obtained H II regions (e.g. the silicon abundance). Moreover, 
a simultaneous characterization of the stellar parameters may address important aspects 
of their evolution, and in particular the role of environmental factors. 

Castro presented the results of a spectroscopic survey in M 33 involving 59 supergiants 
with spectral types between B9 and 09, and a quantitative analysis according to the steps 
described by Castro et al. (2012). A thorough comparison between optical spectra and 
new FASTWIND (Puis et al. 2005) grids resulted in both stellar parameters and chemical 
composition. The parameters derived in conjunction with the evolutionary tracks of Brott 
et al. (2011) hints at the presence of evolved stars with masses in the range 15 and 50 
M0. New routines for deriving the chemical abundances automatically through a process 
optimization showed an oxygen distribution along M 33 that is compatible with previous 
H II region studies (e.g. Rosolowsky & Simon 2008). 

6. Mass loss mechanisms for VMS 

The evolution of Very Massive Stars is presumably dominated by mass loss, which thus 
needs to be understood both qualitatively and quantitatively. Joachim Puis reviewed dif- 
ferent mass-loss mechanisms relevant in this context. Because of the high luminosities, 
only radiation-driven mass loss was considered, and time-dependence, rotation and mag- 
netic fields were not accounted for. 

6.1. Theoretical considerations 
Basic considerations. The equation of motion for the transonic/supersonic regime of 
an expanding wind can be approximated by 

where all quantities have their usual meaning, a is the isothermal sound speed, is 
the total radiative acceleration (lines + continuum), R is the flux- mean opacity per unit 
mass, and T the corresponding Eddington parameter (r e w.r.t. electron scattering only). 
At the sonic point, r s , v = a, and thus = ~g glav implying T(r s ) — 1. To allow for 
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an accelerating wind, g^| s > 0, which then requires ^ | s > and T(r) < 1 below and 
T(r) > 1 above the sonic point, respectively. 

Photon tiring limit. The mechanical luminosity of the wind at 'infinity' is given by 
2 GM\ . ,v 2 v*\ .., / 2GM 



r f Voo WM \ _ 
V 2 R J ~ 



L winii = M^ + ^)=M^ + ^-\ with 



R 



and the maximum mass-loss rate follows from the condition that L w i nd = L* (when the 
star w ould become invisible): M max — 2L lr /{v 00 2 + v csc 2 ). Following lOwocki fc Gavlev 
(|1997l) . M t i r then is the maximum mass- loss rate when the wind just escapes the grav- 
itational potential, with Vqo —¥ 0, and is much larger than typical mass-loss rates from 
line-driven winds, 

M ttI . m . „.„ 32 ^^ * ^ _ „. „ 12 ^r * . 

v csc 2 yr lO 6 !.© i? M yr i? 



Continuum driven winds. To drive a wind by pure continuum acceleration (T = r cont ) 
when the photosphere is sub-Eddington (T(r) < 1 for r < r s ) requires substantial fine- 
tuning to reach and maintain r cont (r) ^ 1 for r ^ r s , and is rather unlikely. Anyhow, in 
such a situation g™^* is almost density- independent, and large mass-loss rates could be 
accelerated, only limited by photon tiring, which needs to be considered in the equation 



of motion (for details, see lOwocki fc Gavlevlll997l ). 



Super-Eddington winds. If, on the other hand, the complete atmosphere is super- 
Eddington, T(r) > 1, continuum driving (mostly due to electron-scattering) might be- 
come possible. When atmospheres approach or exceed the Eddington limit, non-radial 
instabilities arise making them inhomogeneous (clumpy) . Photons on their way out avoid 
regions of enhanced density, and the medium becomes porous. In this case, the photo- 
spheric radiative acceleration decreases compared to an unclumped medium, leading to 
an effective Eddington parameter below unity. In the outer regions, where the clumps 
become optically thin due to expansion , porosity decreases and r^ nt — >■ T(r) > 1. T hus, 



an accelerating wind can be initiated (|Shaviv 2000l 2001al lb). Owocki et al. (2004) ex- 



pressed the effective opacity in terms of a 'porosity length', and showed that associated 
mass-loss rates can become substantial when this length is on the order of the pres- 
sure scale height, amounting to a few percent of the tiring limit. Invoking a power-law 
distributed porosity length, they showed that the 'observed' mass loss from the giant 
outburst of 77 Car might be explained by this (metallicity-independent!) mechanism, and 
that M scales with M oc r 1 / Qp_1 , as long as the exponent of the power-law, a p < 1 and 

r>3. ..4. 



Line-driven winds. The standard theory of line-driven winds (|Castor et al.lll975l and 

later refinements) assumes that the continuum is still optically thin at the sonic point 
(valid for OB-stars, A-supergiants and LBVs in their quiet phase), and that r cont < 1 
everywhere, with r cont — > T c for r ^ r s . The radiative acceleration exerted on a shell of 
mass Am = 47rr 2 /? by one optically thick line can be expressed as 

nnclinc AP L » V dv 



5rad — \ j. a — X 



At Am Aitr 2 pdr ' 

where AP is the transferred momentum, and the term di> arises because of the Doppler- 
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alpha = 0.6 (solid): 0.5 (dotted); 0.4 (dashed) 




Figure 6. M/Mth for a typical O-star wind with v eBC = 600 kms , as a function of r c and a. 

shift within dr. Summing up over all lines using a line-strength distribution function and 
accounting for optical depth effects, the total line acceleration results in 
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with iVo the effective number of driving lines (depending on spectral type and metallicity), 
and < a < 1 related to the slope of the line distribution function. Inserting this 
expression into the equation of motion and neglecting photon tiring, a unique (maximum) 
mass-loss rate can be calculated, 



which dramatically increases for r c — > 1 (in this case, photon-tiring needs to be accounted 
for). On the other hand, the terminal velocity scales with 



N T C \ Va . 
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Figure [5] displays the ratio of M / M t i r as a function of T c , and shows that this ratio is 
strongly sensitive to a. For a = 0.4 (somewhat lower than the typical value for OB-stars), 
the tiring limit would already be reached at T c = 0.8. 



Optically thick winds. The large mass-loss rates from WR-stars, being typically a 
factor of 10 higher compared to OB-star winds at the same luminosity, cannot be ex- 
plained by the standard theory from above. The observed terminal velocities (similar to 
OB-stars) can be reached only when line-overlap effects become efficient. In such dense 
winds, the ionization equilibrium decreases outwards, and photons on their way out can 
interact with lines from different ions, whilst any 'gaps' betwee n lines (as present in 
OB-stars because of an almost frozen-in ionization) are closed fsee lLucv &: Abbottlll993 ; 
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Springmannlll997j) . The initiation of the mass loss, on the other hand, is supposed to rely 
on the condition that the winds are already optically thick at the sonic point, and that 
the (quasi-static) photospheric line acceleration due to the iron peak opa cities around 150 
kK (f o r WNEs) or 40 kK (for WNLs) is sufficient to overcome gravity ( Nugis &: Lamersl 
2002fhlGrafener fc Hamannl (|2005l 120081 ) calculated self-consistent models for WNEs and 
WNLs, and showed that this mechanism actually allows for large M, where the ma- 
j or prerequisite is a high T. Such optically thick wi nds might be pr e sent a lso in VMS 
( Vink et al. 2011 and these proceedings), although IPauldrach et al. ( 20121 ) argue that 
VMS winds might remain optically thin. 



6.2. Monte Carlo mass-loss rates 

Vink discussed three relatively new aspects concerning mass-loss calculations from the 
Monte Carlo method (Abbott & Lucy 1985) - as previously used to predict M for canoni- 
cal OB-type stars (e.g. Vink et al. 2000). The first one concerned the wind dynamics. Until 
2008, the methodology was semi-empirical, as a velocity law was assumed that reached a 
certain empirical v^. Miiller & Vink (2008) suggested a line-force parametrization that 
explicitly depends on radius (rather than the velocity gradient, as in CAK theory), and 
predicted Voo values in reasonable agreement with observations. Muijres et al. (2012) 
tested the Miiller & Vink approach, and as both methods gave similar results, it was 
used in the following. 

Secondly, a new parameter space was probed, i.e. that of the VMS. Vink et al. (2011) 
M predictions show a kink in the M - T relation. For "low" T optically-thin O-star winds, 
the M oc T x relation is shallow, with x ~2, whilst there is a steepening at high T, with x 
~5. At high r the objects show optically thick WR-like winds, with optical depths and 
wind efficiencies above unity. Grafener et al. (2011) provided empirical evidence for such 
a steep exponent (x —5), but there are still issues with the predictions of absolute 
values in this high T range. Critical comparisons between observations and theory are 
underway by Bestenlehner et al. in the context of the VFTS survey (Evans et al. 2011) 

In another study Vink & Grafener (2012) calibrated wind mass-loss rates using an an- 
alytic method to find that the wind efficiency number equals unity right at the transition 
point between optically thin and thick winds: r\ = r = 1. Application of this relation to 
the most massive stars in the Arches cluster suggests there is little room for additional 
mass loss during e.g. LBV eruptions, and current wisdom would suggest that PSN explo- 
sions are unlikely, unless one were to move to lower Z galaxies (e.g. Langer et al. 2007, 
Yoshida & Umeda 2011, Yusof et al. in prep.). 



6.3. Alternative mass loss: eruptions and mass transfer 

Solar metallicity VMS likely evaporate as the result of stellar wind mass loss. However, 
alternative mass loss may also be important, especially for the lower initial mass and sub- 
solar metallicity ranges. Furthermore, we know Eta Car analogs and supernova impostors 
exist in external galaxies (e.g. Van Dyck et al. 2005, Pastorello et al. 2010, Kochanek 
et al. 2012), but quantitative estimates on the integrated amount of such eruptive mass 
loss are hard to obtain as both the eruption frequency, and the amounts of mass lost per 
eruption span a wide range with LBV nebular mass estimates varying from ~0.1 Mq in 
PCygni to ~10M Q in Eta Car (Smith k Owocki 2006). 

The energies required to produce such giant mass eruptions are very high (~ 10 50 erg), 
and their energy source is unknown. Soker (2004) discussed that the energy and angular 
momentum required for Eta Car's great eruption cannot be explained with a single-star 



scenario. 



Very Massive Stars 



137 



There is a growing amount of evidence that the most massive stars are oftentimes 
found in binary systems, and binary evolution with mass loss is pursued by many groups 
around the globe (Vanbeveren 1998; Eldridge et al. 2008). What has become particularly 
clear from recent spectroscopic radial velocity surveys (e.g. Sana et al. 2012) is that there 
is a particularly large number of short-period binaries, which might merge still during 
core H burning, and subsequently evolve as seemingly single stars. For these reasons 
one of the most famous recent phrases in the massive star community has been "binary 
stars might actually be the best single stars" (de Mink et al. 2011). However, also after 
the main sequence, there are still many physical processes involving mass loss through 
Roche lobe overflow (see e.g. Langer 2012) and common envelope evolution (Ivanova et 
al. 2012), which remain as yet ill-understood. 

6.4. Mass-loss diagnostics 

The traditional ways of determining mass-loss rates of (very) massive stars involve re- 
combination lines (such as Ha and He n 4686), as well as radio and sub-mm continuum 
measurements that measure the amount of free- free emission (Wright & Barlow 1975, 
Lamers & Cassinclli 1999). Especially the free-free method may be used in the near 
future with new facilities such as ALMA coming online. The drawback of the above diag- 
nostics is that they depend on an uncertain amount of wind clumping (Puis et al. 2008). 
The unsaturated resonance lines of trace elements, such as P v, located in the far UV, has 
been considered a more accurate tool for mass-loss diagnostics (e.g. Fullerton et al. 2006), 
because their formation depends linearly on density such that inhomogeneitics average 
out. However, it was shown by Oskinova et al. (2007) that the formation of resonance 
lines is also affected by wind clumping, if the line opacity makes the individual clumps 
optically thick. Neglecting this effect may lead to underestimations of the true M (see 
also Sundqvist et al. 2010). 

For these reasons it is important to (i) gain a greater understanding of both the physics 
and the diagnostics of wind clumping throughout the stellar atmosphere and wind, ideally 
as a function of radial distance, and (ii) to develop diagnostics that are not dependent 
on wind clumping. 

6.4.1. X-ray diagnostics for VMS 

Massive stars of most (but not all) spectral types are sources of X-ray emission. In single 
stars, the X-rays most likely originate in the gas heated by the strong shocks resulting 
from the line-driven instability of stellar winds (e.g. Lucy 1980, Owocki et al. 1988). 
Therefore, the properties of X-ray emission are sensitive to the wind driving mechanism. 
Because of the proximity of VMSs to the Eddington limit, the details of wind driving 
and line-driving instability growth may be different from lower mass massive stars. 

The X-ray luminosity of VMS stars is challenging to predict. The X-ray luminosity 
of Galactic OB stars follows the trend Lx oc 10~ 7 Lb o i- While some binary O stars with 
colliding winds have X-ray luminosity significantly higher than 10~ 7 Lb o i) the majority of 
O star binaries follow this correlation as well (Oskinova 2005, Naze 2009). In some cases, 
the binary O stars have X-ray luminosity significantly lower than the expected for a single 
star of similar spectral type. Oskinova produced a diagram showing the dependence of 
the X-ray luminosity in binary O stars on the period. No correlation was seen, and the 
short period binaries can have low X-ray luminosity. 

Oskinova concluded that a low X-ray luminosity cannot serve as a robust argument 
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against a binary nature of an star. In other words, a binary luminous VMS could po- 
tentially have a low Lx- However, the story for the WR-like VMS binaries (with strong 
winds) might be different from the weaker-winded O-star binaries. 



6.4.2. A new wind measurement approach using X-rays from colliding wind binaries 

Sugagawa presented Suzaku observations of the WR binary WR 140, taken at four 
different times around periastron passage in 2009 January. The X-ray spectra changed 
in shape and flux at each phase. As periastron approached, the column density of the 
low-energy absorption increased, indicating that the emission from the wind-wind colli- 
sion plasma was absorbed by the dense WR wind. The luminosity of the dominant hot 
component from the wind-wind collision is not inversely proportional to the (variable) 
distance between the two stars. In the case of the mass-loss ratio Mq/^wr =0.04, Suga- 
gawa could explain this discrepancy if the O-star wind collides with the WR wind before 
it has reached its terminal velocity, leading to a reduction in its wind momentum flux. 
Sugagawa presented these mass-loss rates, which were calculated using the absorptions 
and variations of the spectra (Sugawara et al. 2012, submitted). 

Daminelli showed that the He II 4686 line in Eta Carinae displays two peaks before 
periastron in good correlation with the X-ray intensity. There even is a third peak, in 
close coincidence with periastron, which is anti-correlated with X-rays intensity. This 
may be interpreted as a collapse in the wind-wind collision structure, when most of the 
energy escapes in the extreme UV, which would be possible if the eccentricity is larger 
than e> 0.9. 

What is clear is that the approach presented by Sugagawa using X-ray observations is 
applicable to other massive CWBs with elliptical orbits. In addition, unexpected X-ray 
brightening of very massive CWBs (such as WR 21a) may be helpful for understanding 
VMS mass loss. 



6.4.3. Mass loss at very low metallicity 

Whilst there are still several uncertainties in our empirical knowledge of Solar metal- 
licity mass loss, rates at low Z may provide additional constraints on the driving mech- 
anisms. In this context, Tramper presented the results of a quantitative spectroscopic 
analysis of VLT/X-Shooter observations of six O-type stars in the low- metallicity galax- 
ies IC 1613, WLM and NGC 3109 (Tramper et al, 2011; Herrero et al. 2011, but see also 
Herrero et al. 2012). The obtained stellar and wind parameters can be used to probe the 
mass loss versus metallicity dependence at metallicities below that of the SMC. Tramper 
et al. compared their derived mass-loss rates with the empirical results from Mokiem et 
al. (2007) for the Galaxy, LMC and SMC, and with the theoretical prediction from Vink 
et al. (2001), and argued that the mass-loss rates appear to be higher than expected. 

It is clear that the analysis of a larger sample of stars at sub-SMC metallicities is 
needed to confirm or disprove these results. 



7. VMS Evolution and Fate 

Alexander Heger (with Woosley & Chen) presented the fourth review talk on (very) 
massive star evolutio n. An introduction to massive star evolution can be found in lWooslev et al 



(2002): lHegen(l2012l) . One of the most exciting questions concerning massive stars is how 
they will die, which of them will explode, and how. A range of outcomes is possible in 
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terms of observational signatures, but the key question is how massive the star is at 
the time of death. Therefore, according to Heger understanding the mass-loss rates of 
massive stars is one of the key uncertainties that needs to be resolved. 

Other stumbling blocks along the way involve rotation and binary evolution. For the 
sake of simplicity, Heger only discussed single stars. However, it is clear that some close 
binaries will merge early during their evolution, already during central hydrogen burning. 
These would likely evolve in a fashion similar to that of single stars - though likely rapidly 
rotating. Rotation however is a second parameter that makes the picture more complex. 
So for the sake of simpl icity Heger dealt "just" with single stars and deferred the reader 
to the rotation works of lYoon et al.l ( 20061 120121 ). In some aspects, rapidly rotating single 
stars may evolve in a fashion similar to non-rotating stars of a different (usually higher) 
mass, maybe if mass-loss rates are varied (artificially) similar results can be obtained. 

For the final death, the explosion mechanism and magnitude, knowing the size of the 
(helium) core is the key ingredient. Heger discussed outcomes in terms of initial masses 
for Pop III stars where it was also assumed there is no mass loss by stellar winds, taking 
the key uncertainty out of the equation, and hence allows one to draw a clearer picture 

though obviously not realistic. However, it allows one to understand and formulate the 
possible outcomes. 

Supernovae from lower mass massive stars generally produce neutron stars within the 
'classical' core collapse SN scenario. For higher masses, again really meaning for higher 
mass cores at the end of evolution, one eventually has to deal with core structures that 
no longer allow for an efficient explosion and expect the star to collapse to a black hole. 
Depending on the rotation of the star, a GRB may result. Again, this may require binary 
stars or rapidly rotating single stars, but here the focus is on the dependence on star/core 
m ass for possible outcomes . A recent review on such possible outcomes can also be found 



Wooslev fc Hegerl (|2012h . 



When the mass of the star exceeds some 90 — 100 M Q (~ 42 M He core), pulsational 
pair instability (Puls-PSN) can occur. These are violent nuclear-powered pulsations dur- 
ing the final burning stages, usually powered by oxygen or silicon burning. They may 
produce as much energy as a SN, or more. However, it could also be significantly less. 
Heger discussed this in more detail (see below). After the pulsing is over and the star 
forms an iron core, the fate should be similar to the more massive stars discussed above. 
For higher initial masses, above some 140 M® (65 M Q He core) the first pulse is already 
powerful enough to disrupt the star and a pair instability (PSN) results; for even higher 
initial masses, above some 260 M Q (133 M Q He core) it is expected that the star collapses 
to a black hole. A brief discussion on what could happen in that case may be found in 



Frver et al.l (|200lh . 

One of the striking problems with the PSN theory is that they produce a unique chemi- 
cal signature - at least for the non-rotating primordial stars. Simulation s of the formatio n 
of the first stars suggest that they should have been rather massive (e.g. JAbel et~aTll2002l ). 
hence making PSN, and their nucleosynthesis products should have been incorporated 
into the next generation of stars, but to date, no stars with such nucleosynthesis pattern 
have been found. 



7.1. Pulsational pair instability (Puls-PSN) 

Figure [7] depicts the dynamics of a Puls-PSN simulation in ID: after the pulse there is 
a ring-down phase in which radiative dampening by neutrino losses brings the star back 
into hydrostatic equilibrium. It is noted that the core of the star is much more extended 
after the pulse than it was before (the reader may follow convective boundaries as a 
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Figure 7. Energy loss (purple), energy generation (blue), and convection (green hatching) 
during a pulsational pair instability pulse of a 100 Mq star. The x-axis indicates the time, t, in 
years since the beginning of helium burning, the y-axis shows the radius coordinate from the 
center of the star in solar radii. 

reference). In other words, the energy deposited by the burning during the pulse leads to 
an expansion and cooling of the stars - as we know from textbooks, stars have a negative 
(gravothermal) heat capacity. But in order to get ready for the next pulse, it needs to 
get hot enough, and the energy needs to be lost again. The more energetic, the cooler 
the core of the star will be. But now the end stages of massive stars are dominated by 
neutrino losses, and these are known to depend steeply with temperature. Therefore, the 
cooler the star, the longer it takes to cool. 

In extreme cases of rather powerful explosions, the temperature could drop enough 
such that the core has to cool on the classical Kelvin-Hclmholtz time for radiative losses 
from the surface, as it has become so cool that neutrinos are no longer efficient. This can 
increase the recurrence time by several orders of magnitude. In summary, we generally 
expect that more powerful pulses result in longer delay times, whereas weak pulses result 
in short delay times for the next pulse. The recurrence times can range from hours to 
more than 10,000 yr, and the energies from 0.0001 B to several B (1 B = 1 x 10 51 erg). 

Why is all this relevant? Because we wish to understand whether some giant eruption 
LBVs could be the result of Puls-PSN. In the case of Eta Car we dealt with a fairly weak 
explosion, as some of the H remained bounded to the star, so we would have expected 
the next pulse to happen soon after Eta Car's giant eruption. However, this was not 
observed, or so it appears. Therefore, at this point it seems unlikely that Eta Car's giant 
eruption was a Puls-PSN event. 
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7.2. Stellar Envelope Inflation 

Alternative explanations for LBV outbursts and eruptions have been proposed over the 
years (see Humphreys & Davidson 1994; Vink 2009). Grafener et al. (2012) proposed 
the possibility that envelope inflation near the Eddington limit may play a key role in 
explaining the radius increases during S Dor cycles. The peculiar structure of inflated 
envelopes, with an almost void region beneath a dense shell could mean that many in 
reality compact stars are hidden below inflated envelopes, displaying much lower effective 
temperatures (see also Ishii et al. 1999; Petrovic et al. 2006). 

During the JD, Grafener discussed the inflation effect for WR stars, whose observed 
radii are up to an order of magnitude larger than predicted by theory. Based on a new 
analytical formalism, he described the radial inflation as a function of a dimcnsionlcss 
parameter W, which largely depends on the topology of the Fe-opacity peak, i.e., on 
material properties. For W > 1, an instability limit is found for which the stellar envelope 
becomes gravitationally unbound, i.e. there no longer exists a static solution. Within this 
framework one may also be able to explain the S Doradus-type instabilities for LBVs like 
AG Car (discussed by Groh during the meeting). Moreover, due to the additional effect 
of sub-photosphcric clumping, it may be possible to bring the observed WR radii in 
agreement with theory (see Sander et al. results discussed earlier). 

It should be noted that stellar effective temperatures in the upper HR diagram may be 
strongly affected by the inflation effect. This may have particularly strong effects on the 
evolved massive LBV and WR stars just prior to their final collapse, as the progenitors 
of supernovae (SNe) Ibc, SNell, and long-duration GRBs. 

7.3. 3D Simulations of Thermonuclear Supernovae from VMS 

Ke-Jung Chen (with Heger & Woosley) presented results from numerical simulations of 
the demise of VMS with initial masses between 140 M Q and 250 M Q that can die as 
powerful PSN explosions. Chen et al. used CASTRO, a new multidimensional radiation- 
hydrodynamics code, to study the evolution of PSNe. The 3D simulations start with 
the collapse phase and follow the explosion until the shock breaks out from the stellar 
surface. Unlike the iron-core collapse SNe, PSNe are powered by thermonuclear runaway 
without leaving compact remnants. Much Ni is forged, up to 30 M©, and its decay en- 
ergy powers the PSN luminosity for several months. During the explosion, the emergent 
fluid instabilities cause the mixing of PSN ejecta, and the amount of mixing is related 
to PSN progenitors. The red supergiant progenitors demonstrate strong mixing, altering 
the spectrum and light curves. 

8. Implications 

8.1. Population synthesis models 

The implications for the existence of VMS may be far-reaching, as VMS may dominate 
both the kinetic wind energy input and the ionizing radiation in the Universe. Because 
of higher temperatures at lower metallicities, VMS may be increasingly UV bright. How- 
ever, their higher luminosities might imply higher mass-loss rates and terminal wind 
velocities, which would account for an increased kinetic wind energy. 

Voss presented his recent population models (Voss et al. 2009) which follow both the 
energy - in the form of kinetic wind energy as well as radiation - and the ejection of 
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radio-active isotopes (such as 26 Al) simultaneously. Below 120 Mq, stellar evolution cal- 
culations predict a strong increase in the ejected mass of 26 Al with stellar mass. The 
ejection of 26 Al from stars above this limit has not been studied in detail, but as the 
mass-loss increases dramatically towards higher masses, it is reasonable to expect the 
26 Al to do the same. The short evolutionary timescale of VMSs mean that all the 26 Al is 
ejected 2-4 Myr after the star-forming episode, and if present, the VMSs will dominate 
the signal for 2-5 Myr populations, but due to the decay their signal will be negligible 
for older populations. Comparing the 26 Al signal from 2-5 Myr massive open clusters to 
their 5-10 Myr counterparts is therefore a promising way to probe the evolution of VMSs 
(Voss et al. 2010; 2012). 

8.2. Wolf-Rayet Stars in the Extraordinary Star Cluster NGC 3125-A1 

The powerful radiative and mechanical feedback from very massive stars 100 M ) 
shape the evolution of star-forming galaxies and their environments. Nearby galaxies (< 
10 Mpc) provide excellent laboratories for studying populations of such stars in sufficient 
detail in a variety of astrophysical environments. 

Aida Wofford (with Leitherer and Chandar) studied the massive star populations of 
clusters Al, A2, Bl, and B2 in blue compact dwarf galaxy NGC 3125, which is located 
11.5 Mpc away and has an LMC-like metallicity. It is unclear from past studies if cluster 
Al hosts an extreme population of WR stars. In addition, the WR star populations of the 
other clusters are not well characterized. Wofford et al. obtained HST/STIS 1200-9000 
A spectra of these four clusters, and higher resolution HST/COS 1200-1450 A spectra 
of cluster Al, on which Wofford focused. The STIS spectrum of this cluster shows that 
the equivalent width of Hell A1640 is three times the mean of local starburst galaxies 
(Chandar ct al. 2004) and three times the value of the strongest Lyman Break Galaxy 
(Erb et al. 2010). This suggests that Al must have a large fraction of WR stars relative 
to the number of O stars. Either Al has a top heavy IMF or it contains a few massive 
stars with very strong winds. The COS spectrum of Al shows the strongest O V + 
FeV absorption feature at 1371 A from a starburst in the local universe. The O V line 
originates in the most massive stars and is sensitive to clumping in the stellar wind. The 
analysis of the O V + Fe V using CMFGEN stellar atmosphere models is underway. 

8.3. Nebular He II 4686 emission: an indirect tracer of massive stars at low metallicities 

Shirazi (with Brinchmann) presented a carefully selected sample of 189 star-forming 
galaxies with strong nebular He II 4686 emissions in Sloan Digital Sky Survey Data 
(SDSS) Release 7. They used this sample to investigate the origin of this high ionization 
line in star-forming galaxies where the ionizing continuum almost certainly arises from 
massive stars. The current stellar population models can predict He II 4686 emission only 
for instantaneous bursts of 20% solar metallicity or higher, and only for ages of 4-5 Myr, 
the period when the extreme-ultraviolet continuum is dominated by emission from WR 
stars. 

Shirazi & Brinchmann find however that 83 of the star-forming galaxies (70% at 
oxygen-abundance lower than 8.2) of their sample do not have WR features in their 
spectra despite showing strong nebular He II 4686 emission. Nevertheless, at higher metal- 
licities Hell is always seen with WR features. Shirazi went on to show that the stacked 
spectra of the non-WR Hell emitters do not show WR features either, which suggests 
that the non-detection of WR features in these galaxies is not due to low signal-to-noise 
data, i.e. it is probably real. 

Shirazi proposed that a possible explanation for the discrepancy between the model 
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predictions and the observed data at very low redshifts could be the result of a spatial 
offset between the location of the WR stars and the region where the He II emission arises 
from. Alternatively, as the non-WR Hell emitters appear to be preferentially present in 
younger starbursts, (quasi)-chemically homogeneous stellar evolution could provide a 
possible explanation, as this may lead to higher stellar temperatures, and perhaps result 
in an elevated He II emission even for main-sequence O stars. 

Shirazi & Brinchmann are currently attempting to disentangle these explanations by 
analyzing higher signal-to-noise spectra of a sub-sample of these galaxies that were fol- 
lowed up with the WHT. 

8.4. Very Massive Stars in IZwl8 

I Zw 18 is a blue compact dwarf (BCD) galaxy with the lowest metallicity known (at 1/30 
- 1/50 the solar value), and is therefore thought to be the best local galaxy template to 
galaxies at high redshift. Although IZwl8 is 15-19 Mpc away, Hubble/STIS imagery 
resolves stars in the galaxy. 

Heap showed that UV color-magnitudes diagrams indicate that the most massive stars 
in the northwest cluster of IZwl8 are as massive as 150 M . Heap also showed that 
Hubble/COS far-UV spectra reveal that the mass- loss rates from stars in the NW clus- 
ter must be very low, as only the N V 1240 doublet has a P Cygni profile, and it is very 
weak. The CIV 1549 doublet is resolved with an edge velocity of only ^250 km/s. The 
emission component of the C IV doublet is quite possibly of nebular origin. Most of the 
12 most luminous stars are bluer than ZAMS models, suggesting that the evolution of 
the most massive stars is affected by rapid rotation. A comparison of the observed CMD 
and FUV spectra with new evolutionary models including rotation will yield valuable in- 
formation about chemical enrichment of the ISM, injected energy via ionizing radiation, 
and types of SN explosions in the NW cluster of IZwl8. 



9. Final words 

One of the key science goals for the James Webb Space Telescope (JWST) is going 
to be the identification of the first galaxies with Pop in stars. These objects may have 
been very massive (up to 1000 Mq). At this epoch black hole formation may have been 
more common than at the current time involving solar metallicities. The first couple of 
stellar generations may also have been responsible for the reionization of the Universe: 
an important cosmological epoch that is soon to be probed via the 21cm line with 
instruments such as LOFAR. 

With the E-ELT, individual (very) massive stars may be observed out to the Virgo 
cluster of galaxies in the Hubble flow, at ever larger distances, and for an increasing range 
of metallicities. Basic insights into mass loss at very low Z involve issues such as the self- 
enrichment of metals through mixing (rotation, overshooting, magnetic fields, binarity) 
and mass loss in close proximity to the Eddington limit, which does not necessarily require 
metals, if the winds are continuum driven. We also need to consider the potential mass- 
loss-Z-dependence for the angular momentum evolution and the quest for pinpointing 
the progenitor stars of long-duration GRBs. 

The fate of massive stars is important for our understanding of the chemical enrich- 
ment of the Universe. Whilst the deaths of stars up to 15 M Q now seem to be pretty 
well known; those of VMS up to 300 M Q are as yet a complete mystery. Their fates may 
involve PSN, Puls-PSN, or just normal hydrogen-poor SNe Ibc, either with our without 



144 



Vink et al. 



an accompanying GRB. Key unknown aspects in this quest involve the strength and 
geometry of the progenitor stellar wind. Equatorial winds would remove angular mo- 
mentum during evolution, whilst spherical winds would not. Linear spectropolarimetry 
should become a particularly powerful tool to study the geometry of winds and disks 
during the evolution of (very) massive stars towards explosion. 

Similar arguments can be provided for the evolution of rotation rates during star 
formation. T Tauri and Herbig Ae/Be stars are the optically visible PMS up to ^15 M Q , 
and the formation of such objects is thought to proceed via disks. The latest observations 
and theoretical developments seem to suggest that ever more massive stars may form in 
similar fashions, i.e. via disk accretion, but before such far-reaching conclusions can be 
drawn, large stellar samples over the full IMF are needed. It is particularly important to 
explore the near IR part, in order to diminish the complications by dust extinction. 

GAIA is widely expected to provide key information about the massive star formation 
history of our own Milky Way. This will involve issues of stellar dynamics, cluster forma- 
tion, the role of massive star binarity, and the physical links between massive stars and 
their lower mass (T Tauri) siblings. 

Once we understand the upper IMF and the rotational/multiplicity properties of mas- 
sive stars, we can concentrate on the evolution of the mass loss and rotation properties of 
VMS. Which objects could make PSNe and which ones might produce GRBs? May these 
transient phenomena become star-formation tracers at high redshift? If their luminosity 
functions could be mapped with respect to their redshifts, individual VMS may even 
allow astronomers to constrain galaxy formation models. 

The future of VMS formation and evolution looks bright! 
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